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Abstract. Wc formulate a method of computing invariant 1-forms and structure 
equations of symmetry pseudo-groups of differential equations based on Cartan's 
method of equivalence and the moving coframe method introduced by Fels and Olver. 
Our apparoach does not require a preliminary computation of infinitesimal defining 
systems, their analysis and integration, and uses differentiation and linear algebra 
operations only. Examples of its applications are given. 
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•£} ' 1. Introduction 

o 

The theory of symmetries of differential equations was created by Sophus Lie more 
than a hundred years ago. One of Lie's greatest contributions was the discovery of 
the connection between continuous transformation groups and their infinitesimal ge- 
nerators, which allows one to reduce complicated nonlinear invariance conditions of 
d.e.s under an action of a transformation group to much simple linear conditions of 
infinitesimal invariance - defining equations of symmetry algebra. Lie's method turned 
out to be a powerful tool for studying differential equations, finding their exact solutions, 
conservation laws, etc. 0, [j], f|, [12], [Tj| [19|, |29|, [13|, |3(|. In almost all cases the 
infinitesimal defining equations of the Lie pseudo-groups of symmetries of d.e.s can 
be derived algorithmically. Lie's method requires an integration of (over-determined) 
system of partial differential equations to find a symmetry group of d.e.s explicitly. In 
the last decade methods which do not use an integration but rather extract information 
about structure of symmetry groups directly from their infinitesimal defining systems 
were developed [p3| , [24] , |27| , 28). It was shown how to obtain the dimension of the 



X 



finite Lie group, and in p3| , |24| it was also shown how to find the structure constants 
& k of the symmetry algebra in the finite-dimensional case. In [lG, T^J the method of 



M| was generalized to the case of structurally transitive infinite Lie pseudo-groups. 
Specifically, it was shown how to obtain the Cartan structure equations of the symmetry 
pseudo-group for a system of d.e.s from its infinitesimal defining equations. 

The theory of infinite Lie pseudo-groups was created by Elie Cartan |§] - @. It does 
not use infinitesimal methods and is based on the possibility to characterize an infinite 
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Lie pseudo-group on a manifold M as the set of projections of bundle transformations 
of a principal fiber bundle M x Q — > M, where Q is some Lie group, that preserve 
a collection of 1-forms r l on M x Q. The equations that express the differentials dr l 
through the t 1 and modified Maurer - Cartan forms fi a of the group Q, 

dr { = Al,fi a Ar j + T] k r j At\ 



a.l 



are called Cartan structure equations; they include important information about the 
pseudo-group (see, particularly, OT, Theorem 11.16]). 

In the present paper we apply Cartan's method of equivalence || ^(J and the moving 



coframe method of [1X0, O] to obtain invariant 1-forms of a symmetry pseudo-group 



of d.e.s. Unlike the approach of [Tg, |17[ , the method used here does not require a 
preliminary computation of infinitesimal defining systems and their reduction to the 
involutive form. 

A system 1Z S of differential equations of order s in n independent variables and m 
dependent variables is locally considered to be the subbundle in the bundle J s (£) of s-jets 
of the bundle £ = R"xl m — > R™. A pseudo-group of symmetries Sym(1Z s ) of the system 
1Z S is a subgroup of the pseudo-group of contact transformations of the bundle J s (£) and 
consists of those transformations which preserve the subbundle 1Z S . So the problem of 
finding the group Sym(H s ) is a particular case of the general problem of equivalence of 
embedded submanifolds under an action of a pseudo-group. A powerful and convenient 



moving coframe method for solving this equivalence problem was developed in [ID, [IT 

Some simplifications are possible if we deal with the first order systems of d.e.s. 
By [^2], Theorem 3.3.1.] a system 1Z S is equivalent to the system TZi of the first order, 
which is the subbundle in J 1 (^), where £ = J s_1 (^). The pseudo-group SymiJZi) 
of symmetries of the system 1Z\ is a subgroup of the pseudo- group Cont(J 1 (^)) of 
contact transformations of the bundle J 1 ^)- By Backlund's theorem, ||, [^, Theorem 
4.32], contact transformations on J x {£) are prolongations of point transformations on 
£ . Cartan's method of equivalence allows us to obtain invariant 1-forms which define 
the pseudo-group of contact transformations. Then we can find the invariant 1-forms 
of the pseudo-group Sym{Jt\). To do that, we should make the following steps. First, 
we restrict the invariant 1-forms of the pseudo-group Cont^J 1 ^)) on the subbundle 
Tli and obtain the set of linear dependent 1-forms. Next, we apply the procedure of 
normalization to the appearing conditions of linear dependence. Finally, we apply the 
operations of Cartan's equivalence method to the restrictions on Hi of the structure 
equations of the pseudo-group Cont^J 1 ^)). 
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2. Invariant 1-forms and structure equations of the pseudo-group of 
contact transformations 



According to [521 Theorem 3.3.1.] a compatible system 1Z S of d.e.s of order s is equivalent 
to the system TZi of order 1, which has more dependent variables. So it is possible to 
restrict our attention to the case of s — 1. 

Let 1Zi be a system of partial differential equations of the first order, considered to be 
the subbundle in the bundle of 1-jets of the bundle £ — > X over an n-dimensional 



base manifold X, with g-dimensional fibers. Let 
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...,x n ) denote local coordinates 



of the base X and (u 1 ,^, ...,u g ) denote local coordinates of the fibers of £. Then local 
coordinates of the bundle J 1 ^) are (a? 1 , x n , u 1 , u g ,p\, ■ ■■,pl, ---yPn), an d a 

local section / : X — > £ defined by the equalities u a = f a (x), a G {1, ...,#}, has 
corresponding 1-jet ji(f) : X — > defined by the equalities u a = f a (x), pf = 

ie{l,...,n}. 

A differential form r on J 1 ^) is called a contact form if it is annihilated by all 
1-jets: j'i(/)*t = 0. In local coordinates every contact 1-form is a linear combination of 
the Cartan forms r a = du a — pfdx 1 , a G {1, q} (here and below we use Einstein 

summation convention, so pf dx % = P? dx 1 etc.) 

i=i 

A local diffeomorphism A : J l {£) — > J 1 ^), A : (x,u,p) i— > (sc, is called a 
contact transformation, if for every contact form r, the form A*r is also a contact form; 
in other words, if A*r a = du a — p% doc 1 = £g(x, u,p) t 13 for some functions Q) on J l (£)- 

By Backhand's theorem, ||, ^(| Theorem 4.32], in the case of n > 1 and g > 1 every 
contact transformation A : J 1 ^) — > is a prolongation of a point transformation 

r : £ — ► 5, r : (x, w) t— > (x,u), where the functions pf are defined by the equalities 

„, / dx 1 dx 1 n \ 



du a du a p 
~dx^ + du^ Pj 



Pi 



dxi du 13 



(1) 



To obtain a collection of invariant 1-forms of the pseudo-group of contact transfor- 
mations on J l (£) we apply Cartan's equivalence method |9|, |20|. For this purpose we 
consider the coframe {(r a , dx 1 , dpf) \ a G {l,...,q},i G {1, ...,n}} on A contact 

transformation A acts on this coframe in the following manner: 





( T« \ 




( r« \ 


A* 


dx 1 


= S 


dx 1 




V d P? J 




\ dpf j 



where S : J l {£) — > Q is an analytic function on J x {£) taking values in the Lie group Q 
of non-degenerate block lower triangular matrices of the form 

/ n& \ 





h t ) 



In accordance with Cartan's method of equivalence, we consider the lifted coframe 
on J^S) x Q 

s i = 4e" + 6j^' ) (2) 

where for convenience we use the notations = GLAl, ffa = F?AZ — G%B 3 k d% , 
gfj = Gf k Bj ; (A?) is the inverse matrix of the matrix (ajg), (B{) is the inverse matrix of 
the matrix (ty), so A@ = 8™ and &*■ B 3 k = 5 % k . To find an invariant coframe we use the 
procedure of absorption and normalization of essential torsion coefficients [20, Chapter 
10]. 

Taking exterior differentials of 1-forms 0° and substituting the differentials du@, 
dx 3 , dpj expressed from the equations (2), we obtain 

dB a = (da$ A* + B{ H?° (dj (E2 - fl 6 e - gl E l ) - % E s )) A 7 
+a« B{ H?° S fc AS^- a% B{ E% g\ E k AE l , 

where the functions H? k are defined by the conditions K£ a = 5* <5f . The multipliers 
of E k A E^ and E k A E l are essential torsion coefficients. We normalize them by the 
following choice of the parameters of the Lie group Q: 

Ki = tfBi, (3) 

9% = 9%- (4) 

Then we have 

d& a = §°AQ (3 + E k A E£, (5) 

<S>% = da« A} + c* ffr 0^ - f% E k - c k <% E 3 + eg ££. (6) 

Now the exterior differentials of E l and Ef become 

ds* = a E k + rr 7 a e 7 , (7) 

dZ? = A Sj — J AE« + A£ A 9^ + Og A E 3 , (8) 



where 



n = 4^ fe -4Ef, (9) 

^ = dd 1 + 4^-c k %-c k c^l (10) 

a$ = 4$ + ft, *J + 9% nj - $ - ^ / fc a 7 e e + / fc Q 7 E k 

+c k <% E 3 - c k E«) + (tff + c* E7) + eg f% E7, (11) 

«5 = dflg + & *J + 9% - ./T ; Ef - Ef 

($° - c * /« e' 3 + f^E k + c k g% s E* - c k E£). (12) 
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We note that the conditions (4) imply 

Thus the specifications (3) and (4) of the group parameters of the coframe (2) give the 



lifted coframe 



Q a = a a (du? - dx j ), (14) 



c 



p& + b)<y, (15) 
E? = f% 6" + 9% ^ + a% B\ drf (16) 

on J 1 {£) x 7i, where TC is the subgroup of the group Q defined by the equalities (3) and 
(4). The structure equations (7), (8) do not contain any torsion coefficients, while the 
structure equations (5) contain only constant torsion coefficients. 

The structure equations (5), (7), (8) remain unchanged if we make the following 
change of the modified Maurer - Cartan forms $jg, ^f\, H^, A^, : 

n; ^ n 7 + M| 7 H fe + 7v; 7e e e , 

where if° e , Ljy, M£ 7 , A^ e , f*g 7 , Q^ k , R? jk are arbitrary functions on J 1 (£) x satisfying 
the following symmetry conditions: 

K^e = Ke-yi ^kj = ^>)ki Nye = ^lyi 

-pa pa f^a r^a pa pa pa ( ^'7\ 

r ifi~< — r i-yPi ^i/3k — ^kfiii n-ijk — n ikj — n jik- \ L ' ) 

Their number 

r (1) = \ q 2 (q + 1) + \ n 2 (n + 1) + n 2 q + \ n q (q + 1) + \ n q 2 (q + 1) 
+\ n q 2 (n + 1) + | q n (n + 1) (n + 2) 

is the degree of indeterminancy p(| Definition 11.2] of the lifted coframe O a , H*, 

Using the conditions (|T3|), it is not hard to compute the reduced characters [p0| , 
Definition 11.4] of this coframe: s[ — s' 2 = ... = s q = q + n + nq, s' q+1 = n + nq, 
s ' q +2 = n+ (n-l)q, s' q+2 = n + (n-2)q, ... , s' q+n _ x = n + 2q, s' q+n = n + q, 
s 'q+n+i = s 'q+n+2 = ••• = s' +n+ = 0. It is easy to verify that the Cartan test 



r 



(i) 



s[ + 2 4 + 3 4 + ... + (q + n + nq) s' 



J q+n+nq 

is satisfied, so by definition 11.7 of PUl the lifted coframe (14), (15), (16) is involutive, 
and by theorem 11.16 of PD| , since the last non-zero reduced character s' q+n is equal to 
q + n, the transformations of the invariance pseudo-group of this coframe depend on 
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q + n functions of q + n variables, as it should be. It is easy to verify directly that the 
transformation T : J l {£) x 7i — > J X (S) x 7i satisfies the conditions 

T * e « = e a T*E l = E l , r*Z? = Z« (18) 

if and only if it is projectable on J 1 ^) and its projection A : J 1 ^) — > J l (£), 
A : (x,u,p) t— > (x,u,p), is the prolongation of the transformation r : £ — > £, 
r : (x,u) i— > such that the conditions (Jj) are satisfied. Thus the equalities 

(|T5|) really define the pseudo-group of contact transformations on J 1 (^), when q > 1 
and n > 1. 

Since the forms a , E l , Sf are preserved by the pseudo-group transformations, their 
exterior differentials are preserved also, so T*dQ a = d<d a , T*dE l = dE l , T*dY,f = dEf, 
therefore we have 

T*($° A 9^ + E k A E£) = (T*$£) A 6^ + E k A = $° A 6^ + E k A E£, 

T*(*i a s fc + n 7 a e 7 ) = (T**i) a s fc + (T*n 7 ) a e 7 = % a s fc + n 7 a e 7 , 

T*($° A S 7 - A + A$ A 6^ + fig A E 7 ') 

= T*($^) AS 7 - (T*^ fc ) AEj + (T*A^) A 6^ + (T*fig) A E j 
= $° A S 7 - A Z% + A$ A 9^ + fig A 3' , 

and thus 

T*n = n + ^^ + M^ 7 e 7 , 

r*n 7 = n 7 + m* 7 s fc + /v 7e e e , (19) 

T*A£ = A£ + 7 + Q% k E k + ^ S 7 - M* E£, 
T*fig = fig + + 2 fc - 4 E£ 

for some functions A" 7e , Ljy, M£ 7 , A^ e , P^, Qfp k) Rf jk on J 1 ^) x H satisfying the 
conditions (17). 



3. Symmetries of differential equations 

A suitable method for studying geometrical properties of embedded submanifolds under 
an action of finite-dimensional Lie groups or infinite Lie pseudo-groups was developed 
in [1C, 11]. For its application to the problem of finding symmetries of a system of d.e.s 
1Z\ we restrict the lifted coframe (14), (15), (16) on IZi. That is, we consider the set of 
1-forms 9 a = L*e a , f = i*E\ a? = i*£f, where i : ll x -> J l {£) is the embedding (for 
brevity we identify the map i x id : TL\ x H — > J l {£) x 7i with i : 72-i — > J 1 ^)). The 
1-forms £ l , of are linearly dependent, i.e., there exists a non-trivial set of functions 
U a , Vi, W l a on Tlx x W, such that U a 9 a + V i C i + W l a ^ = 0. 



7 



Setting these functions equal to some constants allows one to express a part of 
parameters ajg, 6* , cjg, gfj of the group H as functions of coordinates of TZ\ and other 
group parameters. Substituting the obtained values of parameters into the modified 
Maurer - Cartan forms (j) a p = ip{ = t***, vr^ = \f p = t*A^, w?. = ^fig- 

makes a part of these forms independent of all differentials of the group parameters. 
Since the transformation Y* changes the forms U.% by the rules (19), in the 

case when the obtained form does not depend on all differentials of the group 
parameters, its coefficients at a J and £ J are lifted invariants of the pseudo-group, and if 
the obtained forms ip l k or ir^ are independent of all differentials of the group parameters, 
their coefficients at crj are lifted invariants also. Normalizing these lifted invariants to be 
constants allows us to express a part of the group parameters as functions of coordinates 
on 1Z\ and other group parameters. If not all group parameters are expressed, we 
should substitute the expressed parameters into the forms 0^, tfj l k , tt*, which depend 
on their differentials, and repeat the process. If the process is completed, but not all 
group parameters are expressed as functions on TZx, we should substitute the modified 
Maurer - Cartan forms <fip, ^1, tv^, Xfp, uj°j, which were reduced during the process of 
normalization, into the reduced structure equations 

d9 a = ^A^ + ^A<, 

dof = A *} - ^ A at + \% A OP + wg A f . 

If the essential torsion coefficients dependent on the group parameters appear, then we 
should normalize them to constants and find some new part of the group parameters, 
which, being substituted into the reduced modified Maurer - Cartan forms, allows us to 
repeat the procedure of normalization. There are two possible results of this process. 
The first one, when the reduced lifted coframe appears to be involutive, outputs the 
desired set of invariant 1-forms which characterize the pseudo- group SymiTZi). In 
the second one, when the coframe is not involutive, we should apply the procedure 
of prolongation pO], Chapter 12]. 

3.1. Example 1: Burgers' equation 

For an application of the above method to finding invariant 1-forms of the symmetry 
group of the Burgers' equation 

u t = u xx + U U X1 

we take the equivalent system of the first order 
u x = v, v x = u t — u v. 
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Denoting x = x 1 , t = x 2 , v = u 1 , u = u 2 , v x = p\, v t = pi, u x = pf, u t = p\, we 
consider this system as a subbundle of the bundle J X {E), £ = R 2 x R 2 -> R 2 , with 
local coordinates {x 1 ,x 2 ,u 1 ,u 2 ,p\,p\,p\,p 2 }, where the embedding t is defined by the 
equalities 

12 12 2 1 

Pi = P2 ~ U U , Pl=U ■ 

The forms 6 a = t*Q a , a G {1,2}, = l* S l , % e {1,2}, are linearly independent, 
whereas the forms = i* E° are linearly dependent. The group parameters ajg, 6* 
must satisfy the conditions det (af) ^ 0, ctet (&*■) ^ 0. Moreover, without loss of 
generality, we can consider that a\ ^ 0, a| 7^ 0, 6} 7^ 0, fr 2 , 7^ 0- Computing the linear 
dependence conditions of forms a? by means of maple, we obtain sequentially the group 
parameters a\, b\, b\, g\ 2 , g\ x , g\ x , j\ 2 , g 2 22 , /| 2 , / 2 2 n as the functions of other group 
parameters and the local coordinates {x 1 , x 2 , u 1 , u 2 ,pl,p 2 } of 1Z\. Particularly, 

2 n 1.2 n i2 b \ a 2 2 (~P2 b 2 + U l U 2 b\ + p^l) <A 

a\ = 0, 6? = 0, b\ = g\ 2 - ~ - 



2 



»1 ' y12 (b}) 3 
a 2 (p1 _ u l u 2 ) l (u l ) 2 a\ - a\p\ - a\p\ - u l (u 2 ) 2 a\ + pla\u 2 + u l u 2 a\ 



9u ~ n.i\2 ' ^11 

/12 = 



{b\) 2 a\ + ^(a}) 2 c 2 &l + uVfai) 2 ^ - u x u 2 a\c\b\a\ - p\{a\) 2 b\c\ + p\a\c\b\a\ 



fli 



{b\fa\ 

u^u 2 {c\a 2 b\a\ - (alYblc 2 ) - p 2 c\a 2 b\a\ + pl{a\?b\c 2 -p\{a\) 2 c 2 b\ + a 2 {b\) 2 

{b\fa\ 



while the expressions for g 22 , f%2 an d /21 are too ^ig to write them out in full here. 

The linear dependence between the forms erf is a\ = a\, a\ = 0. 

The analysis of the modified Maurer - Cartan forms 0^, ip k , 7r* at the obtained 
values of the group parameters gives the following normalizations: 

<t>\ = c\ a\ + (mod 6\ 9 2 , £ 2 , ^2) =► 4 = 0, b\ = % 

^ 2 2 -1^\ = (2c\-c 2 2 )a\ c 2 2 = 2c\; 

ri + ri-^^^claj (mod 9\ 6 2 , £\ e 2 , ^2) 4 = 0; 

2 = - (A\ + ^-^y^ 2 ^ 2 ) e 2 (mod ^ 2 , e, al a!) 

f i alal ~ a\aW + b\(a\) 2 
* hl ~ a 2 

Now the analysis of the structure equations gives step by step the following essential 
torsion coefficients and the corresponding normalizations: 

d6 l = -c\6 2 f\a 2 2 + ... 4 = 0; 
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de l = {{a\ffl 2 - (a 2 ) 2 a 2 + a\a\a\u 2 - {a\) 2 a\b\) 9 2 A ? + ( f 22 + 4 / 2 * 2 A ? + ... 



«2 

,i a2( a 2 Q 2 - a^a^ 2 + ( a \) 2a l b l) f i _ _4 f i 

Jl2- (fl 2 )3 > /22- a 2/21, 



-2 



2 (2 a|a| - a\a 2 u 2 + ^K) 2 ) t 2 x _ aj (a 2 ^ 2 - b\a\) 

(a 2 ) 2 4 A4 +- =► a 2 - 



( (al) 2 (4(a 2 )V-2a 2 ^aX + (a 2 ) 2 ( M 2 ) 2 + (^) 2 (al) 2 ) \ 

(a*) 2 (4 (a 2 ) V - 2 a 2 b\a\u 2 + (a 2 ) V) 2 + (6 2 ) 2 (al) 2 ) 

J21 - /„2\4 



(al) 

, {a\f{b\a\ -a\u 2 ) x 2 x a 2 w 2 _ 

aV 2 = ^ 0A0+... =► & 2 = — , 

do 2 



r 2 = 7^ ^A^... a 2 2 



(a 2 ) 3 s 2 (p 2 2 -v}u 2 yi^ 

jn2 _ 1 n2 , 2 _,_ . 1 _ 1 . 

3al( P i-^ 2 ) 2 /3 y A(7 2 + - <h (p 2 _ u l u 2)2/3' 

mi Z 2 ^ 2wl \ „i , i 3m1 



(p| - u^y/y s " "' y12 (p|- M % 2 ) 2 / 3 ' 

... ' x , 2 (-2 (p 2 ) 2 + 7nVp 2 - 5 (^V) 2 + 2 (^) 3 - 3^p|) \ x 
rf(T 2 = -^22 + 7-2 TTTToTi s A 42 + .. 



(p 2 — M 1 ^ 2 ) 2 

, _ 2 (-2 (p 2 ,) 2 + 7 wVp 2 - 5 {u l u 2 f + 2 (u 1 ) 3 - 3 w 1 p*) 



022 = 



(p 2 — U^U 2 )' 



Thus all the group parameters are expressed as the functions of the local coordinates 
{x 1 , x 2 , u 1 , w 2 ,p 2 ,p 2 } of the equation 1Z±. The result of all normalizations is the invariant 
coframe 

! rf-u 1 — (p 2 — i^u 2 ) dx l — p\ dx 2 
(p 2 — w 1 ^ 2 ) 2 / 3 

A„,2 _ „, 

fl2 



dw 2 — m 1 dx 1 — p 2 dx 2 



^ = (p 2 _ M V)V3 (o^ + u 2 ^ 2 ), 

e 2 = ( P 2 -«V) 2 / 3 ^ 2 , 

x dp 2 — w 2 dp 2 + ((n 2 ) 2 — 2 u 1 ) dw 1 + u l u 2 du 2 

° 2 = {p 2 - M% 2 ) 4 /3 

ii x (p 2 - v}u 2 ) dx 1 + (4 (w 1 ) 3 - 7 (w 1 w 2 ) 2 + 11 u x u 2 p\ - 4 w 1 ^ - 4 (p 2 ) 2 ) dx 2 



+ 



(To = 



{p\ - M % 2 ) 4 /3 

rfp 2 — u 2 du l — u l du 2 — {p\ + v}(u 2 ) 2 — (u 1 ) 2 — u 2 p\) dx 1 



p\ — u l u 2 
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(4 (V) 2 u 2 + {u 2 ) 2 p\ - u 1 ^ 2 ) 3 -u 2 p-3 u x pl) dx 2 



+ 



p\ — u x u 2 



Its structure equations are 

de 1 = i e 1 a + § 1 a a\ + i x a a 2 + £ 2 A a 2 \ 

= A £l + 1 J #2 A £l + 1 ^2 A a 2 + ^2 A | j 



de = e 2 Ae-\eA 



a. 



2' 



d^ 1 = -6* 1 A - 6 I9 l A £ 2 - § I6 2 A C 1 - 2 A a 2 . - 1 A £ 2 

-2 7£ x A o\ + 7£ 2 A a 2 , + f a\ A a 2 , 
da\ = -39 1 A f + 6 2 A - § I0 2 A £ 2 + f 1 A ^ - | J^ 1 A a 2 , 

where the only invariant J is of form 

2 (^ + M V) 2 - (m 1 ) 2 - ?/ 2 p 2 ) 
~ 3 (p| - u% 2 ) 4 / 3 ' 

Taking its exterior differential, we have 

dl = -\0 2 - 2I 2 ? + 2? + - §/a 2 , 

so all differential invariants of the group are functionally expressed as functions of J, 
the rank of the coframe |2(J Proposition 8.18] is equal to 1, and its symmetry group 



is 5-dimensional [20, Theorem 8.22] (as it should be; for the full details of finding 



infinitesimal generators of this group by Lie's method see, e.g., |J0], Chapter 3, § 5].) 

3.2. Example 2: One- dimensional equations of gas dynamics in Lagrange coordinates 

One-dimensional dynamics of polytropic gas in Lagrange coordinates is described by the 
system of d.e.s p6[ 



pt + P 2 u m = 0, 

u t + p m = 0, (20) 
p t + -fppu rn = 0. 

Denoting p = u 1 , u = u 2 , p = u 3 , t = x 1 , m = x 2 and using the above method, we 
obtain the invariant coframe of the symmetry group of the system (20) 

e 1 = -L (du l + (m 1 ) 2 ^ 2 dx 1 - p\ dx 2 ) , 



u 

32 _ 



/ u 1 
7« 3 



(du 2 + pi dx 1 — p\ dx' 



9 3 = — ^ (du 3 + yu 1 u 3 p\ dx 1 - p\ dx 2 ) , 
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II 1 



e = \— 3 dx\ (2i) 

£ 2 = vt-p^dx 1 , 

i i Fr* V\i, (7 - 1) ^) 3 {p 2 ) 2 u 3 - {p\)\u 3 ) 2 - (pp V) 2 2 \ 

^"^VV \ d ^-^ dU 2u^f dX J' 

1 , 7-1 



11 P2 \ ^ 

a 2 = 2 7= i 5 ( d Pl + 7 u 1 p\p\ dx 1 - ^— - u 1 (p 2 ) 2 cte' 
P2 V 7 w w V 2 

(since from considering the physical meaning we have = p > and it 3 = p > 0, 
therefore there is no need to worry about the signs of the expressions under the square 
roots). 

The structure equations of this coframe are 



dO 1 = 


^A^ + ^A^-^Aa 2 , 




d9 2 = 


\ e 1 a e 2 + id 2 a e 3 + h e 2 a e 1 + V ° 2 


A £ 2 + e A a 2 - £ 2 A a 3 , 


de 3 = 


9 1 A £ 2 + 7 2 3 A + ^ A a 3 - £ 2 A <x 2 , 




de = 


±^ A - f# 3 A e 1 - A a 2 , 




de = 






do\ = 


i^^-^e^e-le 1 Ao\-\(2ii- 


- 7 2 + 7) ^ 3 a e 1 



+ A a\ + h e A a\ + 7 (7 - 1 ) ? A o\ - 7 J 2 A a\ + a, 1 A o\ 

da 3 = 6* 1 A + J a 1 A £ 2 - § 6 1 A a 3 - § 9 3 A a 3 + (7 - 1) f 1 A a 2 

-J a ^ A a 3 - I 2 i 2 A a 2 - a 2 A a 3 . 

The invariants 1% and J 2 are defined by the equalities 



7 u 1 - P2 u — P2 u / 7 PI 



1 V U 3 ^{U 1 ) 2 P2 ' ^ V U 1 U 3 p\ 

Their exterior differentials are 

dh = - I j dl + l ( J i " + \o\ - ho\ + l<?l 

so all differential invariants of the symmetry group depend functionally on Ii and I^. 
Thus the coframe (21) has the rank 2, and the symmetry group of the system (20) is 
6-dimensional. In 0, Chapter 3] the explicit form of the infinitesimal generators of this 
group is given. 
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3.3. Example 3: Liouville's equation 

For finding invariant 1-forms and structure equations of the symmetry pseudo-group of 
Liouville's equation 

utx = e u 

we take the equivalent system of the first order 
u t = v, v x = e u . 

Using the notations u = u l , v = u 2 , t = x 1 , x = x 2 and applying the above procedure 
of absorption and normalization, we have a\ = 0, a 2 = 0, while 8 1 , 8 2 , £ 1 , £ 2 , a\ and a\ 
constitute the lifted coframe 





du 1 


— u 2 dx l — p\dx 2 , 




9 2 = 


al( 


du 2 — p\dx l — e" 1 dx 2 ^j 


1 


e = 


(4) 


~ l dx l , 




e = 




11 da; 2 , 




a\ = 


(4 


^e^dpl - (a 2 )~W 


+ a 2 2 gl 2 e ul dx 2 , 


o\ = 


(4 


f (dp\ - vfdx 1 + ((a 2 ,) 


-igly + u 2 p\)dx l ) 



The exterior differentials of these forms are 

dO 1 = -6 2 A + i 1 A a\, 

d0 2 = xi ?\ 2 - 6 1 A^ + ^Aa 2 , 

de = ~Xi^e, (23) 

do\ = X 2 A - Xi A o\ - 6 1 A [a\ + 
rfa 1 2 = x 3 Ae 1 + 2xi A a 2 , 

where 

Xi = {.a 2 2 Y l dal + a\u 2 i x , 

X 2 = dg\ 2 + 2g\ 2 (xi + 1 ) + (a*)-^-" 1 pfe 1 - a\) + w^ 2 , (24) 
Xs = dg 2 u - 3 g 2 uX i + {a 2 2 ) 2 {p{ + (u 2 ) 2 ) (6 2 + f) + 3 a\u 2 a\ + 

W\ and w 2 are free parameters. The structure equations (23) do not contain any torsion 
coefficient depending on the group parameters. The coframe (22) is not involutive, 
because its degree of indeterminancy is 2, whereas the reduced characters are s[ = 3, 
s 2 = ... = s' 6 = 0, so Cartan's test is not satisfied. Therefore we should use the procedure 
of prolongation [^D], Chapter 12]. For this purpose we unite both coframes (22) and (24) 
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into the new base coframe, whereas W\ and w-z turn into the new group parameters. 
Finding exterior differentials of Xii X2 an d X3i we have 



rfX2 = ^iAe 2 -2^Axi-2xiAx2, 

d X 3 = u 2 Ae + 2(9 2 + e) A a\ + 3 xi A X2, 



(25) 



where 



dwi + 3wi(9 +xi 
dw 2 + 4 1^2X2 + 2 



,2\-l -2m 1 



\v\f-9l2) (C 1 + o» - (a 2 2 )-V 



+ e 2 



,2\2 



2\2 



2p 2 



+3a2« 2 x 3 - 



The structure equations (25) admit the change 



for the free parameters 21 and 22- So the degree of indeterminancy of the coframe (22), 
(24) is = 2 again, while the reduced characters now are s[ = 2, s' 2 = ... = s' 9 = 0. 
Cartan's test is therefore satisfied, and the coframe (22), (24) is involutive. Since the 
last non-zero reduced character is s[ = 2, the symmetry pseudo-group transformations 
depend on two arbitrary functions of one variable. This agrees with the result found by 
Liouville In jn| [TjJ the structure equations of this pseudo-group are derived using 
a different method; see also 125 



4. Conclusion 



The approach to computation of symmetry groups used here does not require obtaining 
infinitesimal defining systems, analysis of their involutivity and integration, and includes 
only differentiation and linear algebra operations. So it is algorithmic in principle, 
although the labyrinth of corresponding computations is very intricate. In the future 
it seems that it will be possible to reduce the complexity of computations by means of 
using the canonical contact forms [ETJ on bundles of higher order jets. 
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